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Abstract. We discuss X-ray line formation in dense star winds. A random distribution of wind shocks is assumed 
to emit X-rays that are partially absorbed by cooler wind gas. The cool gas resides in highly compressed fragments 
oriented perpendicular to the radial flow direction. For fully opaque fragments, we find that the blueshifted part 
of X-ray line profiles remains flat-topped even after severe wind attenuation, whereas the red part shows a steep 
decline. These box-type, blueshifted profiles resemble recent Chandra observations of the 03 star £ Pup. For 
partially transparent fragments, the emission lines become similar to those from a homogeneous wind. 
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1. Introduction 

Recent Chandra and XMM-Newton observations of 
O stars show resolved X-ray emission line profiles from 
H- and He-like ions. The line profiles observed exhibit 
a variety of shapes: (1) the stars 8 1 Ori C (Schulz et 
al. 2000) and ( Ori (Waldron & Cassinelli 2001) show 
strongly wind-broadened, symmetric, and nonshifted pro- 
files; (2) <5 Ori A (Miller et al. 2002) and t Sco (Mewe 
et al. 2003) display weakly or nonbroadened, symmetric, 
and nonshifted profiles; (3) ( Pup (Cassinelli et al. 2001; 
Kahn et al. 2001) shows strongly broadened, asymmetric, 
and blueshifted profiles. 

Strong wind-broadening should be a consequence of 
the lines forming in the fast wind. Furthermore, symmet- 
ric and nonshifted profiles indicate line formation in an ex- 
tended and optically thin wind, whereas asymmetric and 
blueshifted profiles indicate line formation in a dense wind, 
where X-rays from the stellar back hemisphere (with re- 
spect to the observer) are more strongly absorbed by cool 
intervening wind gas than X-rays from the front hemi- 
sphere. The most puzzling result of the above observa- 
tions is that a putatively dense wind like C Ori's shows 
line profiles suggestive of an optically thin wind. 

In order to gain a better understanding of the above 
observations, we address X-ray emission line formation in 
a structured stellar wind with thin, aligned absorbing lay- 
ers. We consider a wind consisting of very dense, discrete 
fragments of gas shells that float through essentially empty 
space and are oriented perpendicular to the radial flow 
direction. Such a hydrodynamic structure seems plausi- 
ble for radial wind flows reaching Mach numbers of 100 
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and being subject to strong perturbations. Radiative post- 
shock cooling zones should lead here to highly compressed, 
thin gas sheets lying perpendicular to the flow direction. 

More specifically, a wind with oriented absorber lay- 
ers could result from the combined action of the line de- 
shadowing instability, photospheric turbulence, and the 
Rayleigh- Taylor instability. The line deshadowing insta- 
bility compresses the initially homogeneous gas flow from 
the photosphere into thin, highly overdense shells (Owocki 
et al. 1988). In linear approximation, the line deshadowing 
instability has no lateral component (Rybicki et al. 1990), 
and the resulting flow structure obeys spherical symmetry, 
i.e. is shell-like. However, turbulence- induced velocity fluc- 
tuations at the wind base lead to significant differences in 
the spacing of flow structure (shocks and dense gas shells) 
along neighboring wind rays. Therefore, dense gas patches 
form instead of extended shells. The Rayleigh- Taylor in- 
stability should further fragment these patches. 

We neglect lateral motions (e.g. eddies), and assume 
purely radial flow. The flow shall obey spherical symmetry 
in the statistical sense, i.e., the spatial distribution of frag- 
ments depends on radius, but not on latitude or azimuth. 
In the first part of the present paper, the fragments are 
assumed to be fully opaque. This assumption is dropped 
in the second part. 

X-ray emission originates from shocks associated with 
the dense fragments. Hydrodynamic simulations show that 
the shocks occur almost exclusively on the inner, star- 
ward facing side of the fragments, i.e. are reverse shocks 
(Owocki et al. 1988). To remain general, however, we leave 
it open whether strong forward shocks occur on the outer 
fragment faces, too (Lucy 1982). 
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2. Wind model 

2.1. Assumptions 

Along the lines of this general description, we set up an 
idealized stellar wind with highly compressed gas frag- 
ments oriented perpendicular to the radial flow direction, 
making the following assumptions: 

(1) The gas flow is purely radial, and the radial wind speed 
v(r) = Voo is constant. 

(2) No gas resides outside an outer atmospheric radius 
r = 1. Inside a radius e« 1, the optical depth is infinite 
due to the presence of the central star and dense gas above 
the photosphere. 

(3) Between radii e and 1, all absorbing wind gas resides 
in 2-D absorbers of high density, termed fragments, that 
are oriented perpendicular to the radial flow direction. 

(4) The opening angle of an individual fragment as seen 
from the star is infinitesimal and independent of r. 

(5) The fragments are uniformly random distributed as 
function of the spherical coordinates r, 9, and <f). 

(6) On average, N fragments occur along each radial ray. 

(7) Each fragment is fully opaque. (This assumption will 
be dropped in Sect .01) 

(8) X-rays are emitted from within a sphere of radius x < 
1, with the emissivity scaling as n ~ r ~ 2m ; where typically 
m = 2 (emission ~ density squared). There is no X-ray 
reemission by the absorbing fragments. 

Assumptions (1) and (2) are chosen for calculational 
simplicity. The present paper focuses on qualitative as- 
pects of X-ray line formation, and future work shall deal 
with quantitative modelling of line profiles. Assumption 
(3) is the crucial one, of a wind with oriented absorber 
fragments. Consistent with assumption (1), the frag- 
ments partake in the radial wind expansion (assump- 
tion 4) and have uniform radial distribution (assumption 
5). Assumption (6) will be discussed in Sect. 12.21 below. 
Assumption (7) is chosen to derive a principal limiting 
case, and will be dropped in Sect .01 Assumption (8) ac- 
counts for the fact that the X-ray emitting region may be 
smaller than the region containing absorbing fragments. 
Figure^ shows a sketch of this fragmented wind model. 

2.2. Photon escape 

We discuss assumption (6) and its relation to the total 
radial wind optical depth. The 47r solid angle as seen from 
the star is covered N times by absorbing fragments. One 
may wonder how photons can escape from an atmosphere 
with iV-fold coverage by opaque absorbers. This is due to 
two effects, lateral absorber randomization and sphericity. 
Lateral randomization. One easily proofs that covering an 
arbitrary planar surface of area A with N thin layers of 
opaque absorbing material, then cutting each of the lay- 
ers into small pieces and distributing the latter randomly 
and uniformly within A, an area Ae~ N is left empty of 
absorbers. Therefore, photons crossing through A experi- 
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Fig. 1. Fragmented wind, with N = 10 fragments per ra- 
dial ray. The symbols in brackets designate radii. 
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Fig. 2. Photon escape chanels for N = 3. The dotted 
lines indicate where dense shells got fragmented. The lines 
drawn parallel to the z axis indicate photon escape chanels 
towards the observer. 



ence an optical depth r = N. For such a laterally ran- 
domized absorber distribution, there are still on average 
N absorbers along each ray crossing through A, as is spec- 
ified in assumption (6) above. 

Sphericity. We sharpen for the moment assumption (6) 
to the condition that there are exactly N fragments along 
each radial ray. This would be the case if the flow is purely 
radial, without any lateral motions, as is indicated in 
Fig. |21 Despite the AT-fold coverage, photon escape chanels 
occur like the ones shown in the figure. 

Taking assumption (1) of purely radial flow at face 
value, photons could only escape via sphericity. All central 
rays from the origin r = e, p — > encounter N opaque frag- 
ments, and the optical depth is oo, not N. In the present 
paper we do not interpret (1) in this strict sense, but as- 
sume that lateral fragment randomization has occured. 
This is partly justified by noting that the random radial 
distribution of fragments specified in assumption (5) ap- 
pears as lateral randomization when viewed under an an- 
gle 0. In summary, assumptions (6) and (7) fix the total 
optical depth in radial direction at N. 
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Fig. 3. Basic geometry in the wind. 



3. Opaque fragments 

3.1. Optical depth in fragmented wind 

To calculate the optical depth in radial direction for the 
fragmented wind, consider an arbitrary radius r, infinites- 
imal radial increment 8r, and an area element A of arbi- 
trary shape that is oriented perpendicular to the radial 
direction. A shall be much larger than the average local 
fragment area. The radial optical depth increment is (note, 
r = 1 at the atmospheric rim), 



Sr r = 



5Jk 
A 



NA Sr 

ao~7) 



NSr, 



(1) 



where <5Af, is the area of opaque atomic absorbers in the 
volume A Sr. From Eq. (TjQ), the total radial optical depth 
between e and 1 is N, in agreement with the foregoing 
section. We drop from now on e in nonsingular expressions, 
still using equation signs, i.e. e — * is assumed. 

The observer shall be located at z = oo in standard 
pz coordinates. Introducing the direction cosine p = r ■ z 
with unity vectors r and z, one has v z = pv r and Sz = 
Sr J p, see Fig.|3J The optical depth increment for a photon 
propagating in z direction is, 



St z = N Sz p = N Sr. 



(2) 



The factor p accounts for the reduction in optical depth 
caused by the foreshortening of the fragment area in pho- 
ton propagation direction. Note that this is equivalent to 
an angle-dependent opacity 



X = pN. 



(3) 



Another way to look at Eq. (J2J) is that, for a given pathway 
Sz, N Sr is the average number of fragment hits (assum- 
ing a uniform fragment distribution), which is precisely 
the definition of the optical depth increment. This is elab- 
orated upon in the appendix. The total optical depth, 
r(r, p), from an X-ray emission site specified by coordi- 
nates r and p to the outer boundary of the atmosphere 
along the p ray with p = ryl — p 2 is therefore simply the 
difference of the respective radii r at these two locations, 



r(r, p) 




r) if p > 0, 

2p) if p < 0. 



(4) 



We neglect obscuration due to the central region of radius 
e, which is easy to include in a numerical code but would 
cause lengthy case distinctions in the present equations. 

The conclusion is that r(r, p) is independent of p = 
p(r, p) and therefore of p if p > 0, i.e. for the blue part 
of the line profile forming in gas that moves towards the 
observer. Therefore, the optical depth is the same for all 
photons that originate on the surface of a given front hemi- 
sphere. The reason is that r is the expectation value of the 
number of fragment hits between the emission site and 
the outer boundary of the atmosphere. Since the spatial 
distribution of fragments depends neither on latitude nor 
azimuth, one can arrange them, on average, to spherical 
shells. The same number of shells, N(l — r), is crossed by 
each p ray that starts on the hemisphere with radius r. 
For p < 0, on the other hand, i.e. the red part of the pro- 
file, the number of shells crossed is N{r — p) on the back 
hemisphere, and N(l —p) on the front hemisphere, giving 
in total r(r, p < 0) = N(l + r - 2p) as in Eq. (@J. 

3.2. Optical depth in homogeneous wind 

We compare the above expressions for r with those in 
homogeneous wind. The fragments consist of atomic ab- 
sorbers of cross section cr that have number density no / r 2 
in the homogeneous wind in which the fragments form 
(v = Voo] and we chose a reference radius ro = I). 
The total optical depth of the homogeneous wind in ra- 
dial direction is t = J aon(r)dr. In the limit e — ► 0, 
t = aQUo/e, using equation signs if errors of order e 
are made. The optical depth increment in z direction is 
St z = a n(r)5z = teSzjr 2 . Note that for isotropic ab- 
sorbers, no p factor occurs here in contrast to Eq. (J2J). 
According to Fig.|3 



St z = te — , 
P 



(5) 



where 86 = Szy/l — p 2 jr is the change in latitude over the 
pathway Sz. With p being constant, the z integral is re- 
placed by a 9 integral, giving (e.g. MacFarlane et al. 1991) 



r(r, p) 



U i 
P 



acos^t — asinp), 



(6) 



independent of the sign of p. 

3.3. Flat-topped profile from emitting sphere in an 
optically thin wind 

Consider a spherical shell with outflow speed v and uni- 
form surface distribution of X-ray emitters. There shall be 
no absorption in the volume inside or outside the shell, and 
stellar occultation is neglected. The X-ray intensity mea- 
sured by an observer at infinity is I v = dE u j ' dv du dt dA, 
with the usual meaning of the symbols. Integrating over 
n du gives for the flux, F v ~ dE v /dv. Due to the Doppler 
effect, v = i>Q+pv/vth, where vq is the rest frame frequency 
of the line and Vth the thermal speed. Hence dv ~ dp 
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and dv I dB ~ sinf?. The energy emitted at a given line fre- 
quency v(jx) is dE u ~ ?7 dV^ = 77 dr-rd9-2Trr sin 0. Therefore 
F„ ~ dE u /dO ■ d6/dv ~ r 2 drr/(r). Since any reference to 
/i (or 0) has vanished, the line profile is flat-topped. 



3.4. X-ray line profiles for opaque fragments 



We turn to X-ray line formation in a dense wind with 
absorption. For constant radial wind speed v(r) — v^, 
the surfaces of constant projected velocity pv as seen by 
an observer at infinity are straight cones of opening angle 
acos/Lt. The emergent line profile is given by 



f:: 



drr 2 r]{r)e- T ^\ 



(7) 



with r(r, p) from Eqs. or J5J. Only r depends on p, 
i.e. refers to the explicit cone considered. Since r(r, p) from 
Eq. Q is independent of for > 0, the blue part of the 
line profile from the fragmented wind remains flat-topped, 
despite absorption. This is not true for the homogeneous 
wind. All photons emitted from the surface of a given front 
hemisphere are diminished by the same optical depth (the 
fragment hit expectation value), hence the emission line 
profile stays flat-topped even after severe absorption. Note 
that in a wind with nonconstant outflow speed, the flat- 
topped profiles from individual emitting spheres of radius 
r have different widths [—v(r), +v(r)] in velocity space, 
hence the sum of individual line profiles is no longer a 
flat-topped profile. 

We repeat the physical argument for flat-topped 
blueshifted parts of line profiles despite wind absorption: 
divide the X-ray emitting volume into spherical layers of 
infinitesimal thickness. If the wind would be optically thin, 
the line emission from each of these shells would show a 
flat-topped profile. Consider the front hemisphere towards 
the observer where the blueshifted component of the pro- 
file forms. Because the distribution of absorbing fragments 
obeys spherical symmetry, one can arrange them, on aver- 
age, on spherical shells (with uniformly distributed holes) . 
Therefore, all photons reaching the observer at infinity 
from a given radius r in the front hemisphere encounter, on 
average, the same number of absorber fragments, i.e. ex- 
perience the same optical depth. Hence, the blue part of 
the profile remains flat. 

Figure^] shows line profiles from the fragmented and, 
for comparison, from the homogeneous wind. The lat- 
ter shows the well-known red-blue asymmetry, with a 
rather uniform decline from the blue to the red line wing 
(MacFarlane et al. 1991; Ignace 2001; Owocki & Cohen 
2001). By contrast, the profile for the fragmented wind is 
almost a step function, namely flat-topped over the whole 
blueshifted part, and with steep decline to zero intensity 
on the redshifted part. We suggest that such a profile, after 
convolution with the detector response function, resem- 
bles the blueshifted yet symmetric emission line profiles 
observed from £ Pup, see Fig. 3 in Cassinelli et al. (2001). 
By the same token, its blueshift is too large to be reconcil- 
able with ( Ori's observed emission lines. For Chandra and 
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Fig. 4. X-ray emission line profiles. Full: fragmented 
wind with opaque absorbers; dotted: homogeneous wind. 
Parameters are: N = 10 dense fragments per radial ray; 
inner atmospheric radius e = 0.1; termination radius x = 1 
for X-ray emission; radial optical depth t = 3 of the ho- 
mogeneous wind. The abscissa is in units of the projected, 
normalized wind speed, p = t< z /t>oo. 



XMM-Newton, the line center frequency can be measured 
with an accuracy of roughly 5 mA, which corresponds to 
5v/voo w 0.06 at 1 keV, assuming c/voo — 150. This is 
almost one order of magnitude smaller than the blueshift 
Sv/voc w 0.5 in Fig.0J 

As a conclusion from this section, observed X-ray line 
profiles of £ Pup and £ Ori are qualitatively different. Our 
model of a fragmented stellar wind can only explain the 
former. 



4. Nonopaque fragments 

4.1. Fragment optical depth 

We drop now assumption (7), and account for the finite 
optical depth of fragments. The latter shall form close to 
the photosphere in a rather uniform process, and therefore 
contain roughly the same number of atoms. The radial op- 
tical depth of a fragment is the absorber area per fragment 
area, aon /Nr 2 or te/Nr 2 . In the limit 1/N = dr — > of 
infinitely many, optically thin fragments, the integral over 
dr gives t, as must be. On the other hand, for a finite 
number TV of optically thick fragments, the expression for 
the total radial optical depth becomes, 



dr N 



1 — exp 



-te 



(~Nr 2 



(8) 



This is derived in Eq. i|A.6|) of the appendix, but is also 
intuitively clear: the integrand is the product of the hit ex- 
pectation value of a fragment - which is the optical depth 
increment for opaque fragments, see Eq. Q - times the 
transmission probability for nonopaque fragments. For a 
large number of optically thin fragments, N — > 00, the 
exponential can be Taylor-series expanded to first order, 
giving again if = t. In the opposite limit of opaque frag- 
ments, the exponential term is dropped, giving tf = N as 
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in Sect . 12.21 Turning to moderately optically thick frag- 
ments, we shift the lower integration bound in Eq. JSJ to 
and make a variable substitution y = 1/r. Integrating 
by parts, 



tf/N = 1 - e~ a + v/o7r(l - erf-s/a), 



(9) 



where a — te/N and erf x = 2tt~ 1 / 2 J* e^dt. The average 
optical depth per fragment is t{/N — 0.17 for a = 0.01 
(optically thin fragments), is 0.46 for a = 0.1, and is 0.91 
for a = 1 (tf/N = 1 for opaque fragments). 

Cassinelli & Olson (1979) and Hillier et al. (1993) es- 
timate that the total radial optical depth of the wind of 
C Pup could reach 30 at 1 keV. Time-dependent hydrody- 
namical simulations of O star winds suggest that there are 
of order 10 dense shells. Assuming, then, t = 30, N = 10, 
and e = 0.1, one has te/N = 0.3. Around this value of 
te/N, the line profile shape starts to resemble the profile 
from opaque fragments, see Figs.El and El We must thus 
conclude that assumption (7) of opaque fragments is only 
applicable for rather dense winds. 

4.2. X-ray lines profiles for nonopaque fragments 

According to the above, the optical depth r rp of a single 
fragment located at radius r, with direction cosine p., is 



tc 



Nr 2 fi 



(10) 



in z direction. Here, fi accounts for the longer pathway 
through the fragment when seen under an angle 9. The 
formal integral for the line profile shape becomes 



drr 2 r)(r)e- T ^\ 



r(r,/x) = / dr' N 



1 



6XP (]Vr^) 



(11) 



This cannot be solved analytically, but is easily evaluated 
numerically. 

Figure|5J shows X-ray line profiles at different values of 
the wind optical depth t. From Eq. ©, assumption (7) 
of opaque fragments holds if te/N > 1. In accordance 
with this, the line profile in the bottom panel of Fig.[3 
at te/N — 10, resembles that for fully opaque fragments 
from Fig.0J The same is essentially true for the second 
panel from bottom, at te/N — 1. At smaller t, however, 
the profiles for the fragmented and homogeneous wind be- 
come similar, as is seen in the two top panels of the figure. 

At present it is not clear, which fraction of intrinsic 
X-rays become visible as emergent X-rays. For O stars, 
the latter typically amount to 10~ 7 L, with bolometric 
luminosity L (Sciortino et al. 1990). The wind kinetic 
energy, on the other hand, which is the ultimate reser- 
voir for X-ray energy via thermalization of macroscopic 
flow energy into heat by shocks, comprises for Lv^/lc or 
10 -3 — 10~ 2 L if the wind is close to the single-scattering 
limit, L/c ?« Mvoo, and if < c/100. Therefore, there 
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Fig. 5. Normalized X-ray emission line profiles from winds 
with increasing optical depth (top: optically thin; bottom: 
very optically thick). Full lines show the fragmented wind 
model, dotted lines the homogeneous wind. Parameters 
are N — 10, e = 0.1, x = 1, r\ ~ r~ 4 . The numbers at 
the top right give the total radial optical depth t of the 
homogeneous wind, and the ratios of emergent to intrinsic 
X-ray flux. 

is a large potential margin between intrinsic and emer- 
gent X-rays. The models shown in Fig.[S]all lie within this 
margin. 

4.3. Size of X-ray emitting region 

Hydrodynamic simulations show that shock-compressed, 
absorbing fragments significantly reexpand above w 30 i?* 
due to internal pressure, and the wind gradually becomes 
homogeneous again. On the other hand, X-ray emitting 
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Fig. 6. Similar to Fig.[51 but with smaller X-ray emitting 
volume, x = 0.3. Still TV = 10, e = 0.1. 



sorbing fragments. The line deshadowing instability is ex- 
pected to form strong, X-ray emitting shocks that cool ra- 
diatively, leading to highly compressed fragments of cool, 
absorbing wind gas. If the radiative cooling zone is nar- 
row, one can think of both X-ray emission and absorption 
occuring in the same, spherical segments. The shocks are 
of the braking type, i.e. are reverse shocks on the starward 
facing side of fragments. X-ray emitters on the front hemi- 
sphere, fi > 0, are then hidden behind their natal frag- 
ment, whereas emitters on the back hemisphere are not. 
To treat natal fragment absorption, we have to account for 
the finite optical depth of fragments, since else the whole 
blueshifted part of the line profile would turn black. For 
photons from the back hemisphere, we use r(r, /i) from 
Eq. Ijlljl in the formal integral, however with N dr replaced 
by (N — 1) dr because the natal shell is not hit. For pho- 
tons emitted from the front hemisphere, r(r, /i) is replaced 
by r(r, /i) + r,. M , with r rfl the optical depth of the natal 
fragment from Eq. JTHJ. In Eq. Jill), N dr — > (N — 1) dr 
again, to avoid double counting the optical depth of the 
natal shell. 

FigurcEl shows the resulting line profiles. Since T rfl 
from Eq. (|10fl becomes infinite at fj, = 0, the center of the 
profile from the optically thick, fragmented wind is black. 
This kind of double-peaked X-ray emission line profile is 
not observed. The resolution of the Chandra and XMM- 
Newton gratings is sufficiently high as to resolve such a 
profile shape. 

We are drawn to conclude that either forward shocks 
exist of a strength similar to that of the reverse shocks, 
and prevent the drop to zero intensity; or that the assump- 
tion of spatial proximity of X-ray emitters and absorbers 
in 2-D, oriented fragments is too strong. The former alter- 
native seems unlikely; the latter is rather plausible, given 
that both lateral motions (eddies) and a finite extent of ra- 
diative cooling zones could prevent emitters at fj, > from 
being hidden behind their natal fragment. Therefore, a 
homogeneous distribution of X-ray emitters as in assump- 
tion (8) seems more appropriate than a strict coincidence 
of emitters and absorbers in 2-D, oriented fragments. 



shocks (e.g. as caused by collisions of fast wind clouds 
with fragments) may only reach out to 10 i?*, suggesting 
x = 0.3 in assumption (8). 

FigureEl shows line profiles for this value of x. The ma- 
jor differences when compared to Fig.Elare: (i) the decline 
on the red part of the profile is more gradual now; (ii) the 
blue part of line profiles from the fragmented wind starts 
to become flat-topped at te/N > 0.1 already; (iii) for the 
homogeneous wind at te/N = 1, the ratio of emergent to 
intrinsic X-rays is very small, 4 x 10 -14 . 

4.4. Natal fragment absorption 

We assumed so far that X-rays originate from homoge- 
neous emission within a sphere of radius x. Instead, X- 
ray emission may spatially coincide with the oriented, ab- 



5. Summary 

We consider X-ray emission line formation in a spheri- 
cally symmetric stellar wind with constant radial outflow 
speed. X-rays are assumed to originate from hot gas within 
a sphere of radius x. Cool, absorbing wind gas resides 
in flat, highly ovcrdcnse fragments that are oriented per- 
pendicular to the radial flow direction. The fragments are 
thought to result from radial fragmentation of compressed 
gas shells that form due to the line deshadowing instabil- 
ity. The fragments have infinitesimal opening angle, a uni- 
form spatial distribution, and partake in radial expansion. 
Our results are: 

1. For opaque fragments, the blueshifted part of the 
line profile is flat-topped, even after strong absorption, 
whereas the redshifted part drops steeply to zero. The 
overall profile is therefore essentially flat-topped with a 
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Fig. 7. Similar to Fig-EI but with natal fragment absorp- 
tion included. 



full width of Uqo and a blueshift of Uqo/2. This is in stark 
contrast to the well-known gradual decline from the blue 
to the red line wing for a homogeneous wind. The line pro- 
files from the fragmented wind could explain why observed 
X-ray line profiles from £ Pup are strongly blueshifted, yet 
can be fitted by a Gaussian function. 

2. A blueshift of Woo/2 is too large to be in accord 
with the (essentially) unshifted X-ray line profiles ob- 
served from £ Ori and 9 1 Ori C. We suggest that X-rays 
from these stars do not originate in instability-generated 
wind shocks, but in colliding wind shocks or processes in- 
volving magnetic fields. 

3. The assumption of opaque fragments is applicable if 
te/N > 1 (t the radial optical depth of the homogeneous 
wind between radii e and 1; N the number of fragments 
along radial rays), as could be the case for dense O star 



winds like Q Pup's. On the other hand, with decreasing 
optical depth of nonopaque fragments, line profiles become 
similar to those from a homogeneous wind. 

4. X-ray absorption in natal fragments leads to zero 
flux at line center. This kind of double-peaked profile is 
not observed, which suggests that lateral motions and a 
finite extent of radiative cooling zones prevents a strict 
spatial coincidence of X-ray emitters and absorbers in 2- 
D, oriented fragments. 

Quantitative modelling of X-ray emission line profiles 
of £ Pup accounting for an accelerating wind and instru- 
mental broadening will be subject of a forthcoming paper. 
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Appendix A: A Poisson ring 

We take here a different view at X-ray absorption in a wind 
with oriented absorbers, using Poisson statistics. The two 
objectives of this appendix are (i) to make explicit the 
connection of fragment hits and optical depth, and (ii) 
to obtain an expression for X-ray transmission through 
nonopaque absorbers. We assume N = 1 throughout, and 
generalize to N > 1 at the end of the appendix. In the frag- 
mented wind, consider a ring with radius p centered on the 
z axis at latitude 9. The probability that a photon prop- 
agating at impact parameter p hits an absorber fragment 
along the ring is Sw z — Sr within a pathway Sz. Decrease 
Sz until each fragment that is hit along the ring within 
this Sz fully covers the interval [9, 9 + 59]. Next, divide the 
ring into azimuthal segments of width S<p = S9. Since 
the fragment hit probability Sw z in each of these segments 
is small and the number N$ of segments (or "trials") is 
large, the requirements for a Poisson distribution are met. 
The single parameter of the latter is the expectation value, 
A, here of fragment hits, 



A = Sw z 



2TX[ip 



1 



(A.I) 



using Sr = npS<p/(l — /i 2 ). To give an example, consider 
the ring with p = 1/2 and 9 = 45° (n = l/y/2). On 
average, A = ny2 « 4.44 fragments are hit along this 
ring, in a pathway Sz. The independence of the number of 
hits from N$ is easily understood: if N$ is doubled, 5<p, 69, 
Sz, and Sw z drop by a factor of two, leaving A unchanged. 

The probability ir(k) that k fragments are hit along 
the full ring is, 



k = 0,1,2, 



(A.2) 



Photons at different <f> have different absorption histories, 
but for simplicity we assume that the whole ring is illumi- 
nated by an average intensity (I 1 ). Summing the intensity 



8 



A. Feldmeier et al.: X-ray line emission from a fragmented stellar wind 



1° after absorption within Sz over all segments 5<j), 
J°+/° + ...J^=A^</°) = 

+ tt(1) (P) (N^-l + e^) 

+ tt(2) (I 1 ) - 2 + 2e~ r ^) (A.3) 

+ ... 

+ 7r(AT ) (P) N^e^", 

where r rfl is the fragment optical depth at radius r in di- 
rection acos/z. Inserting l|A.2|) and using the binomial for- 
mula for k <C Nj, (as can be safely assumed for a Poisson 
distribution), it follows that 

M)e -g^-^) . (A.4, 

The sum is the power series expansion of the exponential, 

(1°) = (P) e- x e X ( 1 ' 1 " N * ) = (P) e-M 1 -^^). 

(A.5) 

The case N > 1 is treated by making N copies of the 
atmosphere, and eliminating N — 1 fragments along the 
radial rays in each of the copies. The N copies with one 
remaining fragment per radial ray are arranged like pearls 
on a string, and photons pass subsequently through each 
of them. Hence, for N > 1, 5w z in Eq. 1|A.5|) has to be re- 
placed by N 6w z , giving as optical depth increment along 
Sz, 

5t z = N(l-e- T ^)Sr. (A.6) 

For opaque fragments, this reduces to Eq. of the main 
text. According to Eq. (|A.6() . the partial transparency of 
the fragments leads to a reduction of their effective opacity 
by a factor 1 — exp(— •/>«), i.e. we obtain 

X = fJ.N (l-e- T ^) (A.7) 
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